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We evaluate the tp' to J/?/) ratio in diffractive photoproduction in a hght-cone framework, using 
charmonium wave functions extracted from non-relativistic potential models. Contrary to current 
belief, we find that the best estimate for the ratio is a factor 2 to 5 below the data. The mea- 
sured ratio constrains the distribution of the cc component of the charmonium light-cone wave 
function and indicates that it is more compact than in potential models. We predict that the in- 
elastic photoproduction ratio will be bigger than the elastic one, and will equal that measured in 
hadroproduction. 



I. INTRODUCTION 

Quarkonium production is a hard process in which the 
heavy quarks are produced with limited relative momen- 
tum. This kinematic restriction implies that the standard 
QCD factorization theorem does not apply, ie., there is 
no guarantee that the cross section can be expressed as a 
product of universal parton distributions and a partonic 
subprocess. Quarkonium production is thus sensitive to 
the environment and yields new information about the 
dynamics of hard processes. The abundant data imposes 
severe restrictions on theoretical models (for a discussion 
see Ref. In particular, the recent CDF data |^ on 
charmonium polarization at large p± in pp collisions dis- 
agrees with the color octet model prediction |^ . 

The ratio of cross sections for radially excited states 
is sensitive to the time-scale of the production process. 
If all relevant interactions occur while the quark pair is 
compact, r±_ ^ l/mg, the direct production cross section 
is proportional to |$(0)|^, the square of the quarkonium 
wave function at the origin. Conversely, late interactions 
depend on the shape of the wave function out to its Bohr 
radius, r± ^ 0{l/{asmQ)). The a{'4'')/a{J/ip) ratio is 
thus a very interesting quantity which may give clues on 
the correct production mechanism. 

In forward charmonium hadroproduction the ratio of 
■0' to J/^ cross sections is found to be roughly inde- 
pendent of the kinematics and of the size of the nuclear 



target [^U^]. Its value is moreover consistent with the 
expectation based on the proportionality to |$(0)|^ 
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suggesting that late interactions are unimportant for 
quarkonia produced in hadron fragmentation regions. In 
nuclear fragmentation regions the ratio is measured to be 
smaller than in Eq. (|l|). This is explained by the larger 
break-up cross section of the tp' in late interactions with 
the nuclear comovers. 

Preliminary CDF results [|| show that 52± 12% of the 
T(IS') are directly produced. Given that the T(3S') is 
produced only directly one deduces from the published 
total cross sections that a{3S)/(Jdir{lS) = 0.4 ±0.1, 
consistent with the expectation 0.3 ± 0.05 based on Eq. 
(|l|) for bottomonium. 

Diffractive charmonium electroproduction, 
■>jj{nS)p, is believed to occur via two-gluon exchange 
At large Q^, the size of the quark pair in the virtual pho- 
ton wave function is ~ l/Q ^ l/^asUic), and the cross 
section is predicted to be proportional to |<i>„(0)P. The 
-0' to J/ip ratio at large Q is indeed measured [01 to be 
consistent with |$.0,(O)/$j/^,(O)|2 - 0.5 - 0.6 On 
the other hand, for photoproduction ||ll|] the measured 
value 
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0.15 ±0.034 (g^ = 0) (2) 



is about a factor 3 below the value found at large Q^. 

It was pointed out that due to the moderate value 
of the charm quark mass and a factor from the cou- 
pling of the two gluons, the photoproduction amplitude 
in fact probes the charmonium wave function at a trans- 
verse size r±_ which is comparable to the Bohr radius. 
This reduces the ip' contribution due to the node in its 
wave function. The range of transverse size which is 
probed decreases when the virtuality of the photon 
increases. It is thus possible to measure the shape of the 
charmonium wave function using electroproduction data. 

In Ref. an estimate of the node effect gave a value 
.15 for the ratio (^). 

This value was obtained with harmonic oscillator wave 
functions for the bound states and with the weighted pho- 
ton wave function parametrized as a sum of two Gaussian 
functions ||l^]. We find that the result is very sensitive 
to the parametrization and is thus uncertain to at least 
50%. 

We report here a more quantitative calculation of the 
photoproduction ratio, in a light-cone framework us- 
ing charmonium wave functions obtained from potential 
models. Surprisingly, the calculated ratio (j^) turns out 
to be a factor 2 to 5 below the data. We discuss the 
implications of this for the charmonium wave function. 



II. EVALUATION OF THE RATIO 

The measured ip' / J/rp photoproduction ratio is consis- 
tent with being independent of the photon energy ||ll|,0 . 
At the HI energy the incoming photon fluctuates into 
the cc pair long before the target. We therefore work in 
the high energy regime where the transverse size r± of 
the (S'-wave) cc pair is frozen during its interactions in 
the target and distributed according to the photon wave 
function 
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where x is the light-cone momentum fraction carried by 
the c quark. 

Each exchanged giuon couples to the cc pair with a 
strength proportional to the color dipole length r± (in 
the usual approximation where the gluon wavelengths are 
long compared to r±). The forward scattering amplitude 
is then given by an overlap with the charmonium 
wave function <&, 
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dV_L*^(a;,ri)ri$^(x,r_L)* . (4) 



We consider two models for obtained from the 
non-relativistic Buchmiiller-Tye (BT) and Cornell po- 
tentials, respectively pO[ |. In the non-relativistic limit 



there is a simple relation between the light-cone ampli- 
tude ^^{x,r±) appearing in Eq. (^) and the equal-time 
wave function ^^^{r) given by the Schrodinger equation 
(see, eg., Ref. pq|). In momentum space, 



$^(x,pi) = 
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where 
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and the relative factor is fixed by the normalization con- 
ditions. 
Using 
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{jp]_ + ml f ' 
Eq. (0) reads in momentum space 
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We use the Mathematica program of Lucha and 
Schoberl ||l^ to solve the Schrodinger equation for the 
BT and Cornell potentials. Our results for the cross sec- 
tion ratio (0), 



(9) 



are shown in the first line of Table |. The squared ratio 
of wave functions at the origin shown on the second line 
is the result predicted for highly virtual photons. 
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0.033 


0.070 


i$;,«(o)/$y/«^(o)p 


0.65 


0.64 
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r]_ r± 


0.19 


0.23 



TABLE I. The to Jf-tp ratio calculated with the 
Buchmuller-Tye (BT) and Cornell wave functions for elas- 
tic photoproduction (first line), large electroproduction 
(second line), and inelastic photoproduction (last line). 



III. DISCUSSION 

Eq. (Q) gives a cross section ratio which is an order 
of magnitude smaller than the ratio of wave functions at 
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FIG. 1. The light-cone ^p' wave function $^i{x = ■5,r±) 
(solid line) and the weighted photon wave function r\Ko 
(dashed line) appearing in the integrand of Eq. (^), plotted 
versus r±. The normalizations are arbitrary. 

the origin. This means that the photon wave function, 
weighted by probes the charmonium wave functions 
at relatively large separations. As seen from Fig. 1, r\^^ 
in fact gives similar weights to the tp' wave function in 
the regions below and above the node, leading to a near 
cancellation in the -0' integral. This is the reason for the 
small value of Wj^ and for its sensitivity to the potential. 

In addition, the result for the ■0' overlap integral (^) 
is sensitive to the charmonium wave function in the rela- 
tivistic domain: only 60% of the integral comes from the 
momentum region \p\ < 0.9mc (the corresponding num- 
ber for the J/ip is 95%). This means that the theoreti- 
cal calculation is not reliable for the ip' amplitude, since 
the wave function was obtained from the non-relativistic 
Schrodinger equation. The sensitivity to relativistic mo- 
menta has previously been emphasized in Ref. [ p^ . 

We thus have to conclude that there is no reliable the- 
oretical prediction for the cross section ratio Rjf^ ■ 

The weighted photon wave function r'j_^j probes cc 
pairs with r± < 0.5 fm (c/. Fig. 1). Assuming that 
Eq. (^) is valid in this range implies that the light-cone 
charmonium wave functions we used are incorrect, since 
our result (Table |, first line) is a factor 2 to 5 below the 
data (§). On the other hand, the data may be used to 
constrain the physical wave functions. In particular, 

• The range of the photon wave function narrows 
with Q^, allowing a 'scan' of the wave functions 

• The nuclear target ^-dependence of coherent Jf-ip 
photoproduction is consistent with a'^^ = A°'a'^'^ , 
with a = 4/3 [Q. This implies a small rescattering 
cross section, ie., a small transverse size of the cc 
pairs which contribute to J/ip production. It would 
be important to measure the v4-dependence also in 
coherent ip' production. 

• In the case of inelastic J/^p photoproduction only 



one (Coulomb) gluon is exchanged with the tar- 
get. The overlap integral corresponding to Eq. (^) 
then has a single power of rj_ . This means that the 
charmonium wave function is effectively probed at 
lower values of r±_ than in clastic photoproduction. 
In our calculation using potential model wave func- 
tions the effect of changing the power of r± in the 
overlap integral is quite large, as shown in the last 
line of Table |. 

• The mechanism of charmonium hadroproduction is 
still uncertain, but it is likely that a single gluon is 
exchanged with the target. (This is a feature of all 
present models Just as in inelastic pho- 

toproduction, one gluon exchange implies a single 
power of r± in the overlap integral (^ . This is qual- 
itatively consistent with the measured hadropro- 
duction ratio (|l|), which is larger than Rjl^ but 
smaller than the ratio in large electroproduc- 
tion. 

According to the above arguments, we expect the in- 
elastic photoproduction ratio to be the same as that mea- 
sured in (inelastic) hadroproduction, a{ip')/adir{J/''P) — 
.24. The preliminary data on inelastic ip' photopro- 
duction is still too imprecise for a definite conclusion. 

It is perhaps not so surprising that the potential model 
results for the charmonium wave functions are inaccu- 
rate. In particular, those models take no account of the 
fact that the ip' lies only 44 MeV below DD threshold. 
It follows from the uncertainty principle that the ip' wave 
function contains virtual DD pairs with a lifetime (and 
size) around 4 fm. The photoproduction amplitude, how- 
ever, measures only the |cc) component of the wave func- 
tion. It is quite possible that this component is narrowly 
distributed in r± , while cc pairs at larger separations are 
part of higher Fock states which contain gluons and light 
quarks. 
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